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THE NAVIER- STOKES STRESS PRINCIPLE EOR VISCOUS FLUIDS* 

By Ernst Mohr 



The Navier-Stokes stress principle is checked in the 
light of Maxwell's mechanism of friction and in connection 
herewith the possibility of another theorem is indicated. 



SUMMARY 



The Navier-Stokes stress principle is in general pred- 
icated upon the conception of the plastic body. Hence 
the process is a purely phenomenological one, which Newton 
himself followed with his special theorem for one-dimen- 
sional flows. It remained for Maxwell to discover the 
physical mechanism by which the shear inflow direction is 
developed: According to it, this shear is only "fictitious" 
as it merely represents the substitute for a certain trans- 
port on macroscopic motion quantity, as conditioned by 
Brown's molecular motion and the diffusion, respectively. 
It is clear that this mechanism is not bound to the special 
case of zho one-dimensional flows, but holds for any flow 
as expression of the diffusion, by which a fluid differs 
sharply from a plastic body. If it is remembered, on the 
other hand, that the cause of the stresses on the plastic 
body lies in a certain cohesion of the molecules, it ap- 
pears by no means self evident that this difference in the 
mechanism of friction between fluid and plastic body should 
not prevail in the stress principle as .well, although it 
certainly is desirable in any case, at least subsequently, 
to establish the general theorem in the sense of Maxwell. 
Actually, a different theorem is suggested which, in con- 
trast to that by Navier-Stokes, has the form of an unsym- 
metrical matrix. Without anticipating a final decision 
several reasons are advanced by way of a special flow 
which seem to affirm this new theorem. To make it clear 
that the problem involved here still awaits its -final 
solution, is the real purpose behind the present article. 

*"Uber den Navier-Stokesschen Spannungsansatz fur zahe 
Elussigkeit sstromungen. 11 Luf tf ahrtf orschung, vol. 18,, 
no. 9, Sept. 20, 1941, pp. 327-330. 
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THE PEENOMENOLOCrlCAIi THEOREM OF NAVIEH-STOZES 



The present analysis deals with so-called laminar 
flows of fluids, which also includes gases to the extent 
that their density variations are negligible. No limita- 
tions are suggested by the restriction of the study to 
two-dimensional flows which, however, must alv/ays "be envis- 
aged as being three-dimensional, so that the conventional 
three-dimensional identifications can be maintained. 



Figure 1 shows in the usual notation the stresses 

p , p which in vector notation lead to the stress matrix 
-x -y 

' a x T xy 



(1) 

whence the easily understood differentiation process 



T yx °y 




affords the related force by volume. Then with w = (u,v) 
as the flow velocity of the fluid, the Navier-Stokes the- 
orem states that the stress matrix is proportional to the 
matrix of the deformation speeds 




where p, is a characteristic constant for the fluid, 
termed "viscosity." On computing the force by volume 
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of (2) "by means of the right-hand side of equation (3a), 
it is seen that owing to the continuity equation 



+ $1 = 0 
3x 6y 



(4) 



the second matrix in (3a) at the right, contributes 
nothing. Equation (3) specially postulates the symmetry 
of the stress matrix. As the most important boundary 
condition "the adhesion on fixed walls due to the viscos- 
ity" is cited. 

The same principle (3) is applied to plastic bodies, 
with the difference that (j, is no longer a constant. 

To illustrate: 

1. For the flow indicated in figure 2, termed for 
short, plate flow, the shear inflow direction is according 
to (3): 



T = P. 



du 

ay 



= p. 



d 



(5) 



in which form it had been originally expressed by Hewton. 

2. On the "rigid" rotational flow, (fig. 3, 
speed of rotation), the matrix of the deformation 
disappears identically; hence no shear exists. 



O) = 

spe eds 



MECHANISM OE EH I CT I OH ACCORDING- TO MAXWELL 



The lewton or Havier-Stokes equation is, as seen, 
purely phenomenological : other than that, such as the 
manner in which shear inflow direction of figure 2 comes 
into being, it says nothing. And we are indebted to 
Maxwell who made the important discovery for the case 
that the fluid is a real (ideal) gas, that on the plate 
f low the related shear inflow direction is a result of ' 
the transport on macroscopic momentum, as it arises by 
the diffusion: - since the same number of molecules pass 
from below and above through each unit area of the shaded 
surface y = y 0 in figure 2; while the upper ones involve 
on the average a greater macroscopic velocity, the slower 
layer is subjected per unit time to a certain increment of 
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motion quantity, which, is proportional to (u 3 - u x ) , that 
is, proportional to the speed increase at that point and 
hence according to Newton's fundamental law (according to 
which, the change in motion quantity per unit time corre- 
sponds to the action of a force) equals the effect of a 
shear. 

The calculation involved is briefly reproduced. 

In this Maxwell cone ept ion , the shear (similar to 
pressure) refers to the entrained surface element. The 
result is the flow (fig. 4) , in which the speed increase 
is the same as before. 

With the notation 

c amount of the microscopic, molecule velocity 

n number of molecules per cubic centimeter 

An c number of molecules per cubic centimeter, with a 
speed at interval c c + Ac 

m mass of a molecule 

p gas density; hence m n 

the argument is as follows: 

To define the transport from above, we concentrate 
for the first, on the molecules having a speed at inter- 
val c c + Ac, hereinafter called c molecule, 

for short. From these we eliminate those which arrive 
at an angle in the interval £ . . . . # + A^ toward the 
y axis (fig. 5), and call these c,$ molecules. Such 
molecules per cubic centimeter are 

An c 

-g- sin « A* (6) 

This number, multiplied by the content c cos $ of the 
cylinder of surface length c indicated in figure 5, then 
gives the number AN C of c,# molecules striking the 
unit surface per second: 

An c 

AN C = —5— sin $ bi> c cos $ 
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If such a molecule had its last collision at the 
(oblique) distance Xt it transports at its free passage 
the macroscopic motion quantity 



mu 



= m |\i 0 + X cos * | (8) 



where the value of the pertinent quantities for the dashed 
layer, u Q = 0 , is indicated "by subscript 0. In figure 5, 
where such a molecule is plotted, it should "be noted that 
of c,$ and X only their related projections c», 
and X* are shown in plan view. All the c,$ molecules 
together then transport 



AN, 



m fu 0 + \ =03 4 (p?) \ (9) 



or, using an average value X c for X 

AU C m fu 0 + X c cos t J (10) 

Thus the c, ■& molecules advancing from above give the 
transport t + 

C | V 

T c.a ■ ■^ L£ - sin * A* c cos £ mfu 0 + X- cos $ (§2±) \ (11) 
' 2 I 0 c W/J 

whence the integration along $ from 0 to H gives the 

, 2 
transport from above of all c molecules at 



(12) 



and from below at 



An c J*l 1 . /3u\ 1 



T c = — cm ^- u 0 - - X c Qg-kj. - (13) 



hence leaving an excess t„ : 

c 
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T =T + -T = -i An m c X (14) 
c c c 3 .°"e " c \byj 0 K J 

Integrating with respect to all c molecules and 
substituting the average value cX c for cX c « the looked- 
for transport or shear t is: 

that is, the Hewton formula again (5), and by comparison 
with it for the viscosity the explicit expression 



= 3 P ° Xc 



(16) 



In this calculation Maxwell's distribution law for 
velocities at rest is used, an assumption which is cer- 
tainly permitted so long, as in this case, the macroscopic 
velocity u is small relative to the microscopic c 
(reference 1) . The calculation further indicates that the 
same shear is transported by the hypothetical static unit 
surface of the shaded layer y = y 0 in figure 2. It also 
is not tied to the condition that the velocity profile as 
here has a constant rise: instead of (8) it would then 
afford 



mu 



= p(* 0 



+ X cos 



(If), 



(17) 



the dots denoting negligible terms of lower order. 



PR 03LEM 



In this Maxwellian concept the shear t is therefore 
merely "fictitious," since it simply forms the substitute 
for a certain transport of momentum, which in turn is 
caused by the diffusion, as exemplified in figure 6; dur- 
ing time unit 0 .... 1 a certain number of blank mol- 
ecules depart and a corresponding number of white ones 
arrive. As a corollary the method of continuum mechanics 
is for the present inapplicable because the particle de- 
fined according to the continuum no longer' contains all 
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earlier molecules in the state t = 1J To make the con- 
tinuum mechanics applicable requires first the application 
of ' the correction due to the transition of molecules and 
which exactly results in the fictitious stress. Such an 
analysis is obviously impossible for the corresponding 
flow condition of a plastic body, since in contrast to gas 
it has no diffusion. 



Remembering that the conventional foundation of the 
Havier-Stokes theorem rests on the concept of plastic body 
it seems, in any event, desirable to interpret this general 
equation similar to ITewton's special case in the sense of 
Maxwell. But, on applying the argument associated with 
the greatly schematized (fig. 6) to the arbitrary flow 
(fig. 7), for stress there is obtained in respect to the 
entrained surface element, with d/8n as normal differ- 
ent iat ion : 

Sw 



hence, the unsymmetrical stress matrix 



/Su dv\ 
' 3x dx 



. 9u 
W 



(19) 



which exactly agrees with the first matrix in (3a) at the 
right. According to this there is obtained for this ro- 
tation in contrast to the former a shear inflow direction 

T — |A ^velocity increase j. = jj, cu (20) 

which, similar to the plate flow, every faster layer 
exerts on its slower layer below it. 

Which of the matrices is preferable - the symmetrical 
(3) or the unsymmetrical (19)? Undoubtedly the diffusion 
is the controlling characteristic by which a gas differs 
from a plastic body and it would almost be a miracle if 
this difference had no effect in the stress formula. 



The problem is outlined as follows: 
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Plastic "body 


Gase s 


Cohesion of molecules 


Diffusion rather than cohesion. 


Method of continuum mechan- 
ics directly applicable 


Inapplicable. 


Cause of stresses in the 
cohesion of molecules 


In' the case of the plate flow 
the stresses according to 
Maxwell are a result of the 
transport cf macroscopic mo- 
tion quantity attending the 
dif f us ion . 


The stress matrix is pro- 
portional to the matrix 
of the deformat ion speeds 

* 


? 



So, if it succeeds in extending Maxwell's argument to the 
rotational flow, the question can he answered. This is 
attempted herewith. 



FINDING THE ANSWER 



The procedure is as follows: instead of computing 
the excess t r.ansport or Shear t per unit surface in 
motion as indicated by the shaded area in figure 8, try 
directly for the total shear 2ttt t, which this moving 
surface undergoes. But since this area constantly changes 
into itself by the macroscopic motion, the total shear is 
obviously identical with that total shear to which the 
enclosed shaded fluid mass of figure 8 is subjected by 
the outer fluid; and in this interpretation, however, the 
shaded area may be treated as quiescent J In other words, 
the problem simply involves the excess transport through 
this quiescent area; if such occurs, it corresponds to a 
shearing effect on the shaded fluid mass. Further, vis- 
ualize identical coordinate systems x,y,z as indicated 
in figure 8 (z axis at right angles to plane of drawing) 
plotted along the shaded area in close succession for 
each unit surface. If the transport can be computed for 
such a surface element fixed in space, the task is finished; 



ITAJuA Technical Memorandum if o . 1029 



9 



in the event that the transport is altogether positive, 
the unsymmetr ical equation (19) should he applicable. 

Por the actual calculation of such a surface element, 
that of the plate flow is reverted to, even to the extent 
of using the same notation. figure 9 shows in plan view 
the arrival of a 0,6 molecule from distance A (the pro- 
jected tJ'.X' are shown only): the same conveys the mo- 
tion quantity mu(P) , where u(P) is the speed in point 
P and hence slopes slightly downward. Then the two sub- 
sequent facts permit the reduction to the previous case 
of plate flow: 

1. Put 

u(P = u(P«) + ju(P) - u(P")J, (21) 

for u(P) , since the traces on the average cancel for 
all molecules from the symmetry with respect to the yz 
plane, and it is seen that the calculation can he made 
as if only mu(P") were transported from the particular 
molecule. 

2. Since, further, 

OP" = opi + ... = X cos 6 + ... (22) 

for the lengths 3 OP" and OP 1 (the dots denoting terms 
of the order X , hence are negligible), 

u(P«)= u 0 + X cos 6 (**) + ? .. (23) 

^oy/ 0 

Then, the calculation is obviously the same as before, 

while the speed increase (—) now equals the speed of 
rotation u>. \dy/ 0 

It is shown herewith that the shaded fluid mass per 
unit area in figure 8 is subjected to the shear T = p, tt> 
in flow direction by the outside fluid. The similarity 
of plate flow and rotational flow itself is also plain; 
for both are practically the same in vicinity of the mov- 
ing wall, and in both cases this wall continues to perform 
work to the extent that the molecules advancing at a 
slightly lower macroscopic speed are speeded up again to 
full wall speed and then expelled into the fluid; this 
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work then appears in a corresponding heat. Lastly, the 
equation according the unsymmetr ical matrix (19) is no 
longer in agreement with the so-called Boltzmann axiom 
(reference 2) . 

FURTHER REMARKS 



Having recognized the shear' T as fictitious, the 
same is suspected for the pressure p, and found to be 
a fact. It is proved in two stages: first in quiescent, 
then in moving gas. 

In the case of quiescent gas the pressure p on a 
wall is, as is known, the result of the continual impacts 
of the molecules: a molecule flying with speed c at 
angle d against the wall transmits to it altogether the 
vertical motion quantity 2 c cos and these impacts 

produce, according to the kinetic gas theory, the mean 
pressure 




(24) 



with p density and c 3 the average value of c 3 . A 
hypothetical cut within the gas then leaves the same 
pressure (24) on either side, but in this instance, in 
the Maxwellian sense as transport of molecular or micro- 
scopic speed. In fact, each molecule emerging at angle £ 
with speed c together with the reflective incoming 
molecule furnish an increment of perpendicular quantity 
of motion of the same amount 2 c cos $ as before, hence 
result in the same pressure as on the fixed wall. 

If the gas is in motion our previous concept of shear 
can be supplemented by inclusion of the pressure as follows: 
the transport of the molecules is accomplished in first 
approximation solely by the microscopic speed; but trans- 
ported are 1) the microscopic and 2) the macroscopic speed. 
The second case then affords the previously known shear, 
the first, the pressure. Herein is embodied the expression 
for the free displaceability of the molecules which Euler 
recognized as representative of the pressure and Maxwell 
of the shear. 

Since the new theorem leads to the same motion equa- 
tions, hence also to the same pressure field, the pressure 
resistance in both instances is the same. This also holds 
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for the frictional resistance so far as the body is at 
rest (or in uniform motion, respectively) and the condi- 
tion of adhesion is satisfied. It is sufficient to per- 
ceive the proposition for a surface element of the "body. 
Locating the origin of the coordinates including the x 
axis in it, the y axis perpendicularly outward, the 
adhesion condition gives 



— = 0, ~ = 0 for y = 0 (25) 
dx 3x 

hence, "because of the continuity equation 

^ = 0 for y = 0 (26) 
By 



wherefrom the proposition follows immediately. Differ- 
ences in frictional resistance, therefore, occurs first 
on arbitrarily moving surfaces such, as for instance, in 
the case of our rotational flow. 



Lastly, it is to be observed that, if the new theorem 
proves correct, all transitions, of which the two extreme 
cases have been considered here, prevail between the plas- 
tic body and the fluid. 



Translation by J« Vanier, 
National Advisory Committee 
.for Aeronautics. 
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Figs. 1,2,3,4 
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Figure 1.- Stress pattern. 
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>% Quiescent wall 
Figure 2.- Plate flow. 
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Figure 4.- Plate flow from a reference system, in 
which the dashed layer is at rest. 
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Figs. 5,6,7,8,9 



Kcosd, 




Molecule 



Figure 5.- Definition c-' 

shear as exc -ss 
of transport of macroscopic 
motion quantity (plate flow) 



Figure 6.- Exaggerated 

diagrammatic 
representation of the 
formation of the stresses 
in the case of the plate 
flow. 






Figure 9.- Shear as excess 
of transport of 
macroscopic motion quantity 
(rigid rotational flow). 
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